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We investigate a pseudo-Hermitian lattice system, wliicli consists of a set of isomorpliic pseudo- 
Hermitian clusters coupled together in a Hermitian manner. We show that such non-Hermitian 
systems can act as Hermitian systems. This is made possible by considering the dynamics of a 
state involving an identical eigenmode of each isomorphic cluster. It still holds when multiple eigen- 
modes are involved when additional restriction on the state is imposed. This Hermitian dynamics 
is demonstrated for the case of an exactly solvable 'PT-symmetric ladder system. 

PACS numbers: ll.30.Er, 03.65.-w, 03.75.-b 
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I. INTRODUCTION 

The Hermitian quantum mechanics is a well-developed 
framework because a Hermitian Hamiltonian leads to a 
real spectrum and unitary time evolution, which pre- 
serves the probability normalization. However, a decade 
ago it was observed that a large class of non-Hermitian 
Hamiltonians possess real spectra [l|, Q] and a pseudo- 
Hermitian Hamiltonian connects with its equivalent Her- 
mitian Hamiltonian via a similarity transformation [3|, |j] , 
quantum theory based on non-Hermitian Hamiltonian 
was established |4|-[l3|. In additional, such a framework 
also indicates the preservation of probability normaliza- 
tion if a positive-definite inner product is employed to 
replace the Dirac inner product. Nevertheless, to date 
the interpretation and the measurement in experiment 
of such an inner product are not clear. While the Dirac 
probability (Dirac inner product) can be measured in a 
universal manner. Therefore being an acceptable theory 
of quantum mechanics, the Dirac probability is of central 
importance to most practical physical problems. 

Parity and time-reversal symmetric (PT-symmetric) 
system has attracted much attention due to recent pro- 
gresses on experimental investigations in T^T-symmetric 
optical systems, observation of passive T^T-symmetric 
breaking in passive optical double-well structure 14] and 
observation of spontaneous symmetry breaking together 
with power oscillations in optical coupled system [T^ 
were carried out. Optics offers the rather unique ad- 
vantage in detection the wave function evolution and 
seem to be the most readily applicable [1^, [l3 • In the 
past two decades, general issues of quantum effects in 
quantum systems have proven to be successfully inves- 
tigated in the framework of quantum optical analogy 
based on the fact that paraxial propagation of light in 
optical guided structures is governed by a Schrodinger- 
like equation \l&\. Actually, the intensity observed in 
optical experiment corresponds to the Dirac probabil- 
ity of electric field envelope. It is not guaranteed for 
generic systems that the Dirac probability is preserved 
even when the non-Hermitian Hamiltonian is time inde- 
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pendent. Nevertheless, the violation of the conservation 
of Dirac probability in non-Hermitian system does not 
contradict the Copenhagen interpretation. The implica- 
tions of pseudo-Hermitian system are still under consid- 
eration, peculiar features were exhibited such as double 
refraction, power oscillations, etc. [18|, [l^l following by 
the suggestion of realization of P7"-symmetric structure 
in the realm of optics [20| , while nonreciprocal Bloch os- 
cillation with no classical correspondence was also shown 
in T^T-symmetric complex crystal |2l| . 

We propose a class of non-Hermitian lattice systems in 
this work, the system is composed by a set of isomor- 
phic pseudo-Hermitian clusters, which connected with 
each other in a Hermitian way. We show that in such 
non-Hermitian systems, Hermitian like dynamics can be 
observed, including the property that the time evolu- 
tion is Dirac probability preserving. This is made pos- 
sible by considering the dynamics of a state involving 
the superposition of an identical eigenmode of each iso- 
morphic cluster in general case. In the case of having 
additional orthogonal modes, it still holds when multi- 
ple eigenmodes are involved. This Hermitian dynamics, 
as well as the quasi-Hermitian behavior, are specifically 
demonstrated for the case of an exactly solvable pseudo- 
Hermitian system. 

This paper is organized as follows. In Section [Hi we 
present the model we focus on and its basic properties. 
Section IIIII consists of an exactly solvable example to il- 
lustrate our main idea. Section IIVI is the summary and 
discussion. 



II. HAMILTONIAN AND BASIC PROPERTIES 

A general tight-binding network is constructed topo- 
logically by the sites and the various connections between 
them. As a simplified model, it captures the essential 
features of many discrete systems. Also it is a nice test- 
ing ground for the study of the non-Hermitian quantum 
mechanics due to its analytical and numerical tractabil- 
ity. Much effort has been devoted in recent years to the 
pseudo-Hermitian lattice system |22-34]. The Hamilto- 
nian of the concerned tight-binding network reads as fol- 
lows 
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FIG. 1. Schematic illustration of the concerned networks, (a) 
A lattice consists of three 5-site isomorphic clusters, where the 
different sizes indicate the factor Ac • The dot lines denote the 
similarity-mapping-type Hermitian structure couplings across 
the clusters, (b) A concrete example which is a two-leg lad- 
der. Each rung is a non-Hermitian cluster, (c) Equivalent 
two-band model Eq. H24|) . Here the double-headed arrow de- 
notes the quasi-canonical commutation relations between the 
eigenmodes a = ± for the same cluster. 



region of on-site imaginary potentials |35l] . In the follow- 
ing we will show due to the pseudo-Hermitian clusters 
combined together in a Hermitian way, that there ex- 
ist quantum states obeying Dirac probability preserving, 
even if their profiles cover the imaginary potentials. 

We start with the eigen problem of the Hamiltonian 
Ha- In single-particle invariant subspace, following the 
well established pseudo-Hermitian quantum mechanics 
[ll| - [l3| . we always have 



afla.CT |vac) = Xa^crda^a |vac) 



H„a 



(4) 



and 



Hia^a.a |vac) = Ke^a!^^ |vac) , (5) 



where a G [li -^] and cr G [1, iV^], the operators aa^a and 
tta a have the form 



O'a.a — ^^/(ctO^p Aqct — 7 ^ gi^ttg^l 



(6) 



where 



N 



0=1 a<l3 



Nd 



Ha0 = Kap ^ 0,'l,lO'l3,l + H.C. 



(1) 

(2) 
(3) 



1=1 



which consists of N isomorphic clusters Ha, with each 
cluster has a dimension Nd- The label a denotes the ath 
subgraph of N clusters, and Oa i {cia,i) is the boson or 
fermion creation (annihilation) operator at the ith site 
in the ath cluster. The cluster Ha is defined by the dis- 
tribution of the hopping integrals {XaJij} where Aq is 
real. The set of clusters are isomorphic due the fact that 
they have the same eigenfunctions and spectral struc- 
tures. Note that terms J2a<i3 ^^P i^ self-adjoint since 
Ha0 — HaQ, which describes the Hermitian connection 
between clusters. And such kind of couplings are the type 
of similarity mapping, which is crucial for the conclusion 
of this paper. The total Hamiltonian H is not Hermitian 
when the matrix Jy is not Hermitian. Figure [IJa) shows 
a schematic example. 

In this paper, we consider the case of Ha being pseudo- 
Hermitian, i.e.. Ha is non-Hermitian but has entirely 
real spectrum. Then H is also pseudo-Hermitian in the 
case of real Ac, hence possessing the common exceptional 
point as Ha- In general, a pseudo-Hermitian Hamilto- 
nian does not guarantee the Dirac probability preserving. 
It has been shown that the Dirac norm of an evolved 
wavepacket ceases preserving as long as it touches the 



2^9lafl'<y - 5ll',2_^gtafl<y' - ^aa' 



(7) 



Note that 1/;^}, {s/o-} and {So-} arc independent of a. 
Then the operators a a, a and aa^a are canonical conjugate 
pairs, satisfying 

\a.a,a,a,a' .a'\± — \a.a,a ,0.a' ,a'\± = 0, (9) 

where \-,\± denotes the the commutator and anti- 
commutator. And accordingly, the original Hamiltonian 
can be rewritten as the form 



H — 2_^ Xa^aO-, 



a,(T^Q.cr 



(10) 



a</3,o- 

which has the following subtle features: (i) The ma- 
trix representation of H with respect to the biorthog- 
onal basis {(vac| aa.o-, Aqo- |vac)} is Hermitian, i.e., 
(vac| ac^cr-H"aQ',(T' |vac) = ((vac| ac'^o-'-H"aa,(T |vac)) ; (ii) 
Although it is a non-Hermitian operator, i.e., H ^ H'', 
straightforward algebra shows that 



[aa,a,aa',a']± ^ ^aa' , [Ua^a, la' ,a']± = 0, (11) 

which indicates that although with nonorthogonality of 
the eigenstates as inherent feature of non-Hermitian sys- 
tem, tta.a and a]^, ^, obey quasi-canonical commutation 
relations due to the Hermitian connection structure be- 
tween clusters. This results in a new type of particle 



statistics, that is rarely observed in Hermitian systems, 
thus becomes highly relevant in the presence of non- 
Hermitian terms. 

Considering an arbitrary state in the form 



\^A0)) 



E' 



|vac) , 



(12) 



as the initial state, where ^^ \ca\ =1 and in which only 
the eigenmode a of each cluster is involved. At instant t, 
we have 



l*<T (t)) ^ Yl "^"^ '■^*aa,<T |vac) 



(13) 



In the framework of metric operator theory, H acts as a 
Hermitian system, obeying unitary time evolution in the 
positive-definite inner product [^]. However to date the 
physical meaning of the positive-definite inner product 
is unclear, while the Dirac probability can be measured 
in a universal manner, e.g. Dirac probability of wave 
electric field corresponds to the light intensity in optics 
and is simple to detect in experiment [l6J , therefore Dirac 
norm is of central importance. The aim of this paper 
is to show that contrary to the nonclassical dynamical 
behavior [iSl |2l| , the unitary Dirac probability dynamics 
can also be observed in the pseudo-Hermitian system. 
Actually inserting ^^ ^, ap^^i |vac) (vac| a^.o-' — 1 into 
Eq. dni), we have 



\^<y W) = E ^°'°'I3:<^ l'^^'^) (■^^'^l '^P,'^'^ '"*aa.a |vac) 

a,l3 

= ^ CaUj3aai3,a |vac) , (14) 

a,/3 



where 



Upa = (vac| ap^crS ' aa.a |vac) 



(15) 



is the propagator in the framework of biorthogonal basis 
and satisfies 



/ , UjaU*^ — Sap, 



(16) 



due to the above mentioned feature (i) oi H. Accordingly, 
the Dirac norm has the form 

||$.(t))|^ = (|a>.(t)))^|$.(t)) (17) 

q',/3' a,l3 

a 

where the relation Eq. ([TB)) is applied and the a- 
independent factor A^- can be obtained from 



It follows that although |<I>cr (i)) is not the eigenstate of 
the entire network system, the time evolution is Dirac 
norm-conserving, this is a direct consequence of the 
quasi-canonical commutation relations. The result pre- 
sented here for the evolution of an arbitrary state involv- 
ing an identical isomorphic-cluster-eigenmode provides a 
new way for connecting the pseudo-Hermitian and Her- 
mitian systems. 

It is worth to mention that this probability preserv- 
ing evolution can also occur for a state involving multi- 
ple eigenmodes. This due to the fact that there always 
exist states, which parts belong to different eigenmodes 
are orthogonal in terms of Dirac inner product, hence 
preserve the Dirac probability. For instance, a state in- 
volves two eigenmodes cti and (72, its parts on ci and (T2 
are spatially separated local states with respect to the 
coordinate space a, then the two parts of the state are 
orthogonal in terms of Dirac inner product and the evo- 
lution of such a state is probability preserving since the 
quasi-canonical commutation relations. We will demon- 
strate this point explicitly via the following illustrative 
example. 



III. PSEUDO-HERMITIAN LADDER 

Now we investigate a concrete example to demonstrate 
the application of the previous result. We consider a 
system of a two- leg ladder [Fig. [ijb)], consisting of N 
dimers as pseudo-Hermitian clusters. The Hamiltonian 
reads 



N 



N 



Hhudd — 2_^ ^c + Z-^ Ha,a + 1, (19) 

a— 1 a—1 

Ha = -J{aa^iaa,2 + H.c.) -t- i7 {na,i - na,2) , (20) 

2 

Ha,a+i = -nY^{a\^^aa+ij + H.c), (21) 

1=1 

where Ua/ — o-a gCLa/ is the particle number operator 
and the operators obey the periodic boundary condition 



with ^ = 1, 2. K (J) is the hopping 



(vac| a'l,^aj3^^ |vac) = A^Sa^- 



(18) 



"^N+i.e — "i, 

integral along legs (rungs) and 7 denotes the norm of 
the imaginary on-site potential. Note that the ladder 
is a T'T'-symmetric Hamiltonian, where V is the parity 
and T denotes time-reversal. The simple structure of 
this model makes it an ideal testing ground for a more 
profound understanding of the Hermitian dynamics in a 
pseudo-Hermitian system. Taking the transformations 

aa,. = -=L= (e-o/'al, - ae^'^/^al J , (22) 
V2cos0 ^ ^ 

Qa.o- ^ —r ^ I e ' aa,\-oe ' aa,2] , (23) 
V 2 cose' ^ ^ 

where (a G [1, A^] , cr = ±), which are obtained from the 
solution of the dimer (a general solution of A^-dimension 



cluster is shown in Ref. [3g].), the ladder Hamiltonian 
can be written as 

N 

Hhadd = 2^ {~Kaa,aa.a + l,cr ~ Kaa+i^cyO-a^a (24) 
a— l,cr— ± 



+aAaa,craa,cr)- 



which is illustrated in Fig. [TJc), here A — \/ J"^ — 7^ and 
sin0 = 7/ J, B G [0,7r/2]. The biorthogonal structure of 
the solution for a dimer admits the following canonical 
commutation relations Eq. ([5]) and 



[al^„,aa',cr\^ = [aa,(T,aL,crJ± =sec9Saa', (25) 

[al_^,aa\cy]± = [aa,a,al^-a]± = icT tSiuOSaa' ■ (26) 

Obviously, Hamiltonian Eq. (j24p represents a two-band 
model, which has an interesting feature comparing to a 
Hermitian two-band model: although there are no inter- 
band transitions, the two bands are not independent. It 
is due to the pseudo-Hermiticity of the clusters, which 
allows [aa,a,aa-cr]± = but [aa,cr,a'l-a]± ¥" 0. This 
characteristic will be further demonstrated through the 
following quasi-canonical commutation relations Eq. ((33)) 
and the time evolution for various Gaussian wavepackets. 
Figure [Ijc) schematically illustrates such an equivalent 
two-band structure. Nevertheless, Hamiltonian Eq. ([M)) 
can be diagonalized as a Hermitian one, i.e., we have 



k,a=± 

ek,± — ~2k cos fc ± a, 
by using the linear transformations 



Ofc.CT — 



Ofe., 



e ■'a,- 



N 
N 






(27) 
(28) 

(29) 
(30) 



where k = 2mr/N, n e [1, A^]. The linearity of the trans- 
formations allows 



[ak,cr,ak',a']j. — Skk'Scrcr', 



0. 



(31) 
(32) 



However, when dealing with the Dirac inner product, the 
quasi-canonical commutation relations 

["L<Tiafc',o-]± = [ak,a,a\,^^]± ^secdSkk', (33a) 



[al,^^,ak',a]± = [ak,a,al, ]± 



icr tan0(5fcfc',(33b) 



will be taken into account. Such quasi-canonical com- 
mutation relations reflect the subtle features of the sys- 
tem: when dealing with different fc, ak,a and a],, ^, act as 
canonical conjugate pairs and the system displays Her- 
mitian behavior. 



We can gain some insight regarding the role of the 
quasi-canonical statistics. We will see shortly that such 
a model displays the similar dynamics as a Hermitian 
ladder. We start our investigation from the quantum 
dynamics of various initial wavepackets. In the situation 
of a Hermitian ladder, any two wavepackets propagate 
independently and the total probability is preserving. 

Considering an arbitrary state involving both upper 
and lower bands 



l*(0)) = Yl h,aak,a\y&c), 



(34) 



k,<T=± 



"^ith X]fe,a=± l/*::crl = 1, we have 

ll'J' {t))? = E l/fc'-l' (^^^1 [4,a:afe,'.]± Ivac) (35) 

k.(T 

+ E A* -^/fc.-e-'2.At ^^^^1 [al_^^-akA± |vac) 

k,(7 

= sec0 + itan0E'^/fcVaA.-e-'"'"(*/^°\ 



where T-q = tt/A denotes the period of the oscillation. 
The first term gives the contribution from single band, 
while the second term captures the influence of the non- 
Hermiticity. For vanishing 6 we recover the unitary evo- 
lution in Hermitian system. Evidently, ||$ {t))\ — sec 9 
for a state with f^ _„fk.a — 0; which involves only a sin- 
gle mode. Note, however, that mathematically speaking 
the time dependent terms can vanish even in the case 
of fl -afk,cr 7^ 0, e.g. additional orthogonality of the 
wavepacket with multiple eigenmodes. To demonstrate 
this, we study the evolution of initial wavepackets of the 
form 



|^(7Va,A^b,, 



n 



^.,0))^^ 



(36) 



g {k-4,AY/{2p^)^-t(k-<t,A)NA^^^ 



+ ^-(k-4,Bf /{^p") ^-l{k-4,B) 



Ne 



a-kr 



which is the superposition of wavepackets A and B, where 



n^2j:^e- 



-{k-^AY/p' 



2Efce"(''"'^«^ /P ■ The time 



evolution of wavepacket is a powerful tool for under- 
standing the dynamical property of Hermitian quantum 
systems [33] ■ Recently, The propagation of wavepacket 
in discrete systems has been utilized as flying qubit for 
quantum state transfer |38l - l42| . In the Hermitian case an 
initially Gaussian state stays Gaussian as it propagates 
for a long time, especially for the case of |0a b| = t^/2 

El- 

For a sufficient broad wavepacket {p ^ 1), we have 
V, « pN/^/n. Equation ([55| can also be expressed in the 
coordinate space spanned by {a]^ ^ |vac), a]^ ^ I'^^c)} as 



0.02 


0.02 


0.02 



;^ 0.04 

=■ 0.02 

% 
,^ 0.02F 




0.02 


0.02 



t=o.oo : 




0.45 : 




0.90 : 


, 


^ 1.35; 




1.80 : 




2.25 : 




2.70 : 





0.04 



/\ 


t=0.00 


1 ' ' 


0.04 



/\ 


0.29 


1 ' ' 


0.04 
? 


XX 


0.58 




?^ 0.04 
S 


0.75 


, 


_- 0.04 



-'-'"V 


0.84 


1 , 1 


0.04 



./x 


0.91 


1 ' ' 


0.04 



/\ 


1.00 


1 



100 200 300 400 

J 



100 200 300 400 

J 



FIG. 2. (Color online) The Dirac probabilities Pi {j, t) (black dashed line) and P2 (j, t) (red solid line) of a particle, initially 
located in the state I* [Na, Nb,cJ>a, 0s, 0)) for a system with N = 400, 7 = 0.05, J = 0.10, k. = 1.00, and p = 0.05. We obtain 
6 — 7r/6 and the time t is in units of Tu ~ 36.276 k~^ . We plot the Eq. (|40|) for two cases with (a) (J)a ~ ~4>b = 7i"/2, Na = 100, 
Nb ~ 300 and (b) (pA = (pB = 7!"/2, Na ~ Nb ~ 150. The shapes of all the curves are in agreement with the analysis in the 
text. 



\-^{NA,NB,(l)A,(f>B,0)) 



-p^U-NAf/2 i4,AJ 



4y^cos( 



{e-'^Ki 



(37) 



.-^e/2i 



a=l 

+ ^-p'U-Nn)V2^^^sJ (^e-^'^/^l.+e^'/^l,^] |vac) , 

which involves both eigenmodes (a = ±) and actually 
composed of four wavepackets with centers at A'^^th and 
TVsth sites of the legs 1 and 2, and with the velocities (f>A 
and ^B , respectively. To investigate the dynamics of the 
Dirac norm, substituting 



fi 



1 



fc,+ 



g-(fc-0A)V(2p')g-»(fc-0A)A'A^ ^gg^^ 



■^'- = vu' 



-(k-4,BY /(2p^)^~l(k~4,B)NB 



(38b) 



into Eq. ([35]). we have 



I* {NA,NB,<i)AAB. i))r = sec0 + 



Pa-4>bY /{ip^) ^-p^iNB-NAY /i 



(39) 
sin {2TTt/T-£) — (Pab) tan 9. 



where ipAB = {Na + Nb) {4>a ~ 0b) /2. We note that 
if the two wavepackets of Eq. (1551) are well separate in 
A; or a space initially (wavepackets orthogonal in k or 
a space), the weighted exponential factor becomes zero, 
then the probability is always conserved in the evolu- 
tion even they meet each other in the coordinate space 



a. This indicates that for states having additional or- 
thogonal modes, Hermitian like behavior still holds even 
multiple eigenmodes are involved. 

To show more detailed propagation behavior, we study 
the profile of Pg (j, t) (£=1,2), where 



Pi{3,t) = |(vac|aj,^|*(7V^,7VB>A>i3,t))|^ 



(40) 



It is a convenient way to investigate the dynamical prop- 
erties from two typical cases: (a) (pA = —4'b = 7r/2, 
\Na-Nb\ > 2Vh[2/p and (b) cj^a ^ (J^b = tt/2, 
Na = Nb- In case (a), the situation corresponds to two 
counter-propagating wavepackets, with the evolved wave 
function 



\^iNA,NB,7T/2,-7T/2,t)) 



n- 



Vn 



X > je " — e 

k 



I*' {Na + 2Kt, Nb - 2nt, 7r/2, -7r/2, 0)) , 



(41) 
I vac) 



where the approximation of Taylor expansions for cos k 
around ±7r/2 are used for two wavepackets and |^') rep- 
resents the superposition of two wavepackets as state 
l^t) but with different overall phases. It shows that 
the evolved state is still the independent nonspreading 
wavepackets. Similarly, the evolved wave function for 
case (b) has the form 



|vI/(iV^,7V^,V2,vr/2,0) = ^E[' 



-(fc-^/2)V(2p^) 



It has more clear profile in the coordinate space £, i.e. 

(43) 



\^{NA,NA.^/2,Ti/2,t)) 

N 



J2 9iit)J2^ 



-p'lj-{NA+2Kt)]^/2 ij^/2 i 



«=1,2 



J = l 



where 



9iit) 



/TTCOSI 



aj^^ |vac) 



cos (vri/To - 61/2) , ^ = 1 
i sin {nt/TD + 9/2) ,£^2 



(44) 

Obviously, it represents two breathing shape-invariant 
wavepackets propagating along two legs of the ladder 
with the breathing period Td. Furthermore, the Dirac 
norm P/ = Y.j Pe U, t) (^ = 1,2) and Pj^ = Pf + P^ can 
be obtained as the form 



PI 

Pi 
p^ 



cos'' {TTt/Tn-e/2)/ cose. 
sin^ {irt/Tu + 9/2) / cos 0, 
sec 6* + tan 6* sin (27ri/TD) . 



(45) 
(46) 
(47) 



As mentioned in the introduction, the profile of the 
evolved wave function Pg (j, t) can be observed in experi- 
ment. In practice, the quantum-optical analogy has been 
employed to visualize the dynamics in the non-Hermitian 
system [18l - l20l |. In this context, the light intensity corre- 
sponds to Pi {j,t) (for a review, see [la]) and the profile 
corresponds to the light intensity distribution along its 
propagation direction. 

It follows that a manifestation of the non-Hermitian 
nature of i?Ladd is represented by the relative phase 6 
between the breathing oscillations of the two legs, which 
also leads to the time-dependent Dirac probability. The 
profiles of the evolved wave functions and the Dirac 
norms are plotted in Figs. [2] and [3l We can see that 
in case (a) the evolved wavepackets propagate indepen- 
dently and the Dirac norms are preserving. It indicates 
that although the Hamiltonian is non-Hermitian, due to 
the quasi-canonical commutation relations which is a di- 
rect consequence of the Hermitian connection structure 
between clusters, it acts as a Hermitian ladder for some 
initial state. In contrast, the dynamics of case (b) differs 
drastically from the Hermitian case and the Dirac norm 
is no long preserved. Further, the phase difference be- 
tween the breathing oscillations on the two legs can also 
be observed in case (b). 




t (units of T ) 



FIG. 3. (Color online) The Dirac norms Pf (t), Pi (t) (blue 
dotted line) and P^ (^) (black dashed line) for the case of 
(j>A = -<t)B = vr/2 [as in Fig. Oa)]. The Dirac norms Pf (t) 
(red triangle), PI (t) (green square) and P^ (i) (black circle) 
for the case of (J)a = 4>b = 7r/2 [as in Fig. [2Ib)]. All the 
parameters are the same as in Fig. [2l The phase difference 
9 = 7r/6 and also the quasi-canonical commutation relations 
sec6 f» 1.155 are indicated . The shapes of all the curves are 
in agreement with the analysis in the text. 



consist of a set of isomorphic pseudo-Hermitian clusters 
combined in a Hermitian manner, that Hermitian like 
dynamics could be observed in such non-Hermitian sys- 
tems, including the property that the time evolution is 
Dirac probability preserving. As an application, we in- 
vestigate a concrete network, a T'T-symmetric ladder, 
composed of many pseudo-Hermitian dimers. It is shown 
that it acts as a Hermitian system in the following sense: 
besides the reality of the spectrum and probability pre- 
serving, the propagation of certain wavepackets exhibit 
the same behavior as that in a Hermitian ladder. Our 
finding indicates that the reality of the spectrum as well 
as the Dirac probability preserving dynamics can occur 
in a system that violating the axiom of Hermiticity. This 
will pave the way for the development of descriptions of 
quantum system and provide a topic of considerable in- 
terest in a wide range of subjects. 
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